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ABSTRACT 

We present the first detailed spatio-kinematical analysis and modelling of the planetary 
nebula Abell 41, which is known to contain the well-studied close-binary system MT 
Ser. This object represents an important test case in the study of the evolution of 
planetary nebulae with binary central stars as current evolutionary theories predict 
that the binary plane should be aligned perpendicular to the symmetry axis of the 
nebula. 

Deep narrowband imaging in the light of [Nii]6584 A, [Oiii]5007 A and 
[Sn]6717-|-6731 A, obtained using ACAM on the Wilham Herschel Telescope, has 
been used to investigate the ionisation structure of Abell 41. Longslit observations of 
the Ha and [Nil] 6584 A emission were obtained using the Manchester Echelle Spec- 
trometer on the 2.1-m San Pedro Martir Telescope. These spectra, combined with the 
narrowband imagery, were used to develop a spatio-kinematical model of [Nil] 6584 A 
emission from Abell 41. The best fitting model reveals Abell 41 to have a waisted, 
bipolar structure with an expansion velocity of ~40 km s~^ at the waist. The sym- 
metry axis of the model nebula is within 5° of perpendicular to the orbital plane of 
the central binary system. This provides strong evidence that the close-binary system, 
MT Ser, has directly affected the shaping of its nebula, Abell 41. 

Although the theoretical link between bipolar planetary nebulae and binary cen- 
tral stars is long established, this nebula is only the second to have this link, between 
nebular symmetry axis and binary plane, proved observationally. 

Key words: planetary nebulae: individual: Abell 41, PN G009.6+10.5 - circumstellar 
matter - stars: mass-loss - stars: winds, outflows 



1 INTRODUCTION 

It has long been believed that planetary nebulae (PNe) 
result from an early dense stellar wind, originating from 
the progenitor asymptotic giant branch (AGB) star, being 
swept into a thin shell by a later fast wind from the emerg- 



* Based on observations made with the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 
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ing white dwarf fWDV lKahn fc WestI (Il985l) l ater e xtended 
this model, first developed by IXwok et all l| 19781 ). show- 
ing that an equatorially-enhanced AGB wind would lead 
to the formation of an aspherical nebula. This model went 
on to become known as the generalised in teracting stellar 
winds (GISW) model (see e.g., review by iBalick fc FrankI 
though GISW is not always consistent with ob- 
servation of_ex£aJision_velocity versus ionization structure 
(see e.g..l Meaburn et al.ll2005l V The mechanism causing this 
anisotropic mass-loss is not yet clear. The most widely ac- 
cepted source of this anisotropy is for the central star of the 
planetary nebula (CSPN) to interact with a binary partner 
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Figure 1. Deep ACAM-WHT images of A 41 in the light of (a) [Nii]6584 A, (b) [Oin]5007 A, (c) [Sii]6717+6731 A and (d) [Nii]6584 A 
divided by [Oni]5007 A, where black indicates a large ratio. Display scales have been chosen so as to highlight structural features referred 
to in the text. 



l|de Marcdl2009l l. To produce the most extreme bipolar PN 
stru ctures (e.g. that of MyCn 18, the Etched Hourglass Neb- 
ula, iBrxce^TaLlIlHg^), it is believed that the CSPN would 
have to form a common env elope (CE) with its binary com- 
panion (|Bond fc Lividll990l ). For a CE to form, one compo- 
nent of the binary must accrete mass from its Roche lobe- 
filling partner more rapidly than it can thermally adjust 
to the additional material, thus also filling its Roche lobe. 
After this point, any further mass lost by the donor will 
go on to form a CE surrounding both stars, and it is this 
CE which goes on to form the high equatorial density re- 
quired by the GISW model |Nordhaus & Blackman 200i). 
This mechanism is seen as the most likely method by which 
very close binaries (for example, cataclysmic variables) are 



formed (jCrauer fc Bondlll983l ). because as the common en- 
velope is ejected by transfer of angular momentum, the bi- 
nary spirals in, dramatically shortening the period of the 
system. 

Abell 41 (PN G009.6-)-10.5, a = 17'^29'"02.03' 
S = -15°13'04.4" J2000) discovered by [Abell fc GoldreichI 
was classified by iBond fc Livid (Il990l) as elliptical 
under the classification scheme of Balick " i[l987l ). However, 
deeper Ha-|- [Nil] 6584 A imagery reveals "that the nebular 
morphology exhibits an 'H' shape with the addition of 
fainter material forming a continuous loop" jPoUacco fc Belli 
IT997i ). 

Photometric analysis of the CSPN, MT Ser, revealed 
it to be a close binary, showing minima at regular intervals 



Abell 41- shaping of a PN by a binary central star 3 



-15:13:00 




03.5 03.0 02.5 17:29:02.0 01.5 01.0 00.5 
Right ascension 

Figure 2. A deep ACAM-WHT image of A 41 in the light of 
[Nll]6584 A showing the positions of the 5 MES-SPM longslits. 



of 2''43™ l|Grauer fc Bondlll983l ). iBruch et all (|200ll ) con- 
firmed the binary nature of MT Ser but were unable to accu- 
rately determine its orbital parameters because they found 
two different models which fit the observed data, (a) The bi- 
nary consists of a hot sub-dwarf and a less evolved secondary, 
in which case the period is 2'' 43™ and the variations are due 
to a refiection effect (inclination, i = 42.52° ± 1.73°) B (b) 
The binary consists of two evolved, hot sub-dwarfs with a 
period of 5''26™' where the variability results from partial 
eclipses and ellipsoidal variations {i = 65.7° ±0.9°). They de- 
termined the optimum parameters for each model, but con- 
cluded that only radial velocity observations would be able 
to distinguish between the two. Subsequent observation and 
modelling by Shiman skii et al.] (2008) confirmed the pres- 
ence of two sub-dwarf components, but gave no i ndepe ndent 
confirmation of the orbital inclination. lOglozal l|2002l ') pre- 
sented photometric and radial velocity observations along 
with modelling of MT Ser, independently determining an 
orbital inclinati on of 72° ± 15°, whi ch is consistent with the 
second model of iBruch et al] (|200ll . i = 65.7° ± 0.9°). How- 
ever, the values determined for the rest of the system param- 
eters differ vastly, indicating that any agreement should be 
treated wit h caution. Add itionally, the stellar temperatures 
derived by lOglozal l|2002l ) are inconsistent with the obser- 
vational detection of tw o hot sub-dwarf binary components 
Shimanskii et al. I2OO8I ). Of the two mod els of lBruch et al.l 



,200ll ) a nd the mode l of lOglozal (|2002l ). only the second 
model of lBruch etHI (|200ll . ^ = 65.7° ± 0.9°) is consistent 
with photometric observations and the detection of two hot 
sub-dwarf central stars, indicating that this is the most re- 
liable modelling of the binary CSPN system. 

In this paper we present longslit spectroscopy of A 41, 
from which, combined with narrowband images, we derive 
a spatio-kinematical model of the nebula, with the aim of 
understanding the relationship between the nebula and MT 
Ser. 



^ Here, the incUnation, i, is defined sueh that for i = 90° the 
orbital plane would be in the line of sight (i.e. eclipsing). 



2 OBSERVATIONS AND DATA REDUCTION 

Narrowband images of A 41 were acquired using ACAM 
combined with the 4.2-m William Herschel Telescope on 
2009 August 4 ([Nil] 6584 A and [Qui] 5007 A) and 2009 
August 29 ([Sll]6717+6731 A). The seeing for both sets of 
observations did not exceed 0.9". ACAM was employed in 
standard imaging mode without binning resulting in a pixel 
scale of 0.25" pixer\ and using the [NII]6584/21, Taurus 
5009/15 and [SII]6727/48 filters (ING fihers #85, #108 and 
#86). Three 15 minute exposures were taken in each filter, 
the data were bias-corrected, flat-fielded and cleaned of cos- 
mic rays using starlink software, starlink software was 
also used to remove the background lunar contamination, 
arising due to the close proximity of A 41 to the Moon, 
from the [Sll]6717-H6731 A images. The resulting images 
were then co-added and are shown in Figure [1] 

Spatially resolved, longslit emission line spectra of A 41 
were obtained with the second Manchester Echelle Spec- 
trometer combined with the 2.1-m San P edro Martir Tele- 
scope (MES-SPM, iMeaburn etHI |2003| ). MES-SPM was 
used in its primary spectral mode with a narrow- band 90 A 
filter to isolate the Hq and [N 11] 6548 and 6584 A emission 
lines of the 87th echelle order. Observations took place on 
two separate runs in 2004 June and 2007 June, both with 
the same instrument and set-up, using a SITe3 CCD with 
1024 X 1024 24 fim square pixels {= 0.31" pixel" ^). AU 
integrations were of 1800 seconds. 

Binning of 2 x 2 was adopted for all the spectral ob- 
servations, resulting in 512 pixels in the spatial direction 
(= 0.62" pixel"^) and 512 pixels in the spectral direction 
(= 4.79 kms'^ pixel"^). The slit used was 30- mm long 
{= 5') and 150 ^m wide {= 2.0" and 15 kms'^). 

Data reduction was performed using starlink software. 
The spectra were bias-corrected and cleaned of cosmic rays. 
The spectra were then wavelength calibrated against a ThAr 
emission-lamp. Finally, the data were rescaled to a linear 
velocity scale (relative to the rest wavelength of [Nll]6584 A, 
taken to be 6583.45 A) and corrected to show heliocentric 
velocity, Vh^i- 

In total, five integrations were obtained, two with slits 
running North-South, one with a position angle (PA) of 
61.8°and a further two with PA equal to 151.8°. The slit 
positions are shown in Figure [5] and the spectra are shown 
in Figure [3] 



3 ANALYSIS 

3.1 lonisation structure 

Figure [T] shows A 41 in the light of [Nll]6584 A (a), 
[Olll]5007 A (b) and [Sll]6717-F6731 A (c). In aU three 
emission lines the nebula displays rough ly the structure re- 
marked upon by iPoUacco fc Belli l|l997l ), an ellipse where 
the western and eastern edges (parallel to the major axis 
of the apparent ellipse) appear much brighter than the rest 
of the nebula. The northern side of the nebula also appears 
brighter than the southern side at all three wavelengths. 
There are, however, differences between its appearance in 
the three emission lines; for example the ring-like feature at 
the c entre of the nebula (t he horizontal in the H analogy 
from IPoUacco fc "Bdil |l993)), is much more well defined in 
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Figure 3. Comparison of observed [Nll]6584 A PV arrays and synthetic equivalents: (a) observed and synthetic PV arrays of Slit 1, (b) 
observed and synthetic PV arrays of Slit 2, (c) observed and synthetic PV arrays of Slit 3, (d) observed and synthetic PV arrays of Slit 
4 and (e) observed and synthetic PV arrays of Slit 5. The slit positions are as shown in Figure [5] and the velocity axis on all plots is 
heliocentric velocity, Vj,ei. 
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Figure 4. PV arrays of Slit 3 (as shown in Figure[2]l in the light of 
(a) H« with white contours marking the 95th, 97th and 99th per- 
centiles and (b) [Nil] 6584 A convolved to the same thermal width 
as that of Ha (assuming = IQ-* K, the original [Nn]6584 A PV 
array is shown in Figure [3jc)). 



[Nii]6584 A and [Sn]6717+6731 A than in [Oiii]5007 A. From 
the imaging alone it is unclear whether this ring is actually 
a material ring or a projection effect, where the ring is ac- 
tually where the two lobes overlap in the line-of-sight. It is 
also unclear whether the lobes themselves are open-ended 
or closed, as although their edges appear brighter (indica- 
tive of a closed shell, due to more nebular material in the 
line-of-sight), this effect is also seen in open-ended nebulae 
(such as MyCn 18: Sahai et al. 1999) 

The nebula appears comparatively brighter centrally 
in the light of [Olll]5007 A compared to [Nll]6584 A and 
[Sll]6717-|-6731 A as shown by the ratio image in in Fig- 
ure Hid). Upon further inspection, this central bright re- 
gion in [Olll]5007 A appears to be an internal bipolar struc- 
ture, although it is not clear whether this is just a line-of- 
sight brightness effect or, in fact, a separate internal nebu- 
lar shell such as thos e seen in other nebulae (e.g. NGC 7009: 
Sabbadin eral]|2004l and Hen 2-104: ISantander-Garcfa et alj 
20081 ). We note that this feature can also be fo und in the 
Ha-|- [Nil] image shown in [Miszalski et al. I (|2009l . originally 
shown in iPollacco fc Bellll 19971 ). 

The nebular position-velocity (PV) profiles also appear 
different in the light of [Nll]6584 A compared to Ha, even 
once the different thermal widths of the two lines have been 
taken into account (see Figure [4]). This ionisation stratifica- 
tion is a common feature among planetary nebulae, where 
the [Nn]6584 A emission is seen to trace the outside of the 
nebular shell wherea s the Hq is relativ ely evenly distributed 
throughout (see e.g.. lGurzadvanlll997l ). 

The nebular PV arrays also show the same brightness 
variations across the nebula as the imaging, for example Slit 
1 [FigurefSfa)] shows a bright partial velocity ellipse from the 
northern part of the nebula but much fainter emission from 
the southern part. Slit 3 [Figure [3{c)], which lies approxi- 
mately parallel to the major axis of the apparent nebular 
ellipse, shows that the northern end of the nebula is red- 



shifted with respect to the south indicating that the nebu- 
lar inclination is such that the northern side of the nebula 
is pointed away from the observer and the south towards. 
The Slit 3 echellogram also shows two separate emission re- 
gions from the near and far sides of the nebular shell (i.e. no 
closed velocity ellipse), suggesting that the nebula is open 
at both ends of this symmetry axis. This is supported by 
the emission profile from, the slit perpendicular to Slit 3, 
Slit 5 [Figure (H^e)], which shows a partial velocity ellipse. 
The ellipse is not closed as a result of both the open ended 
and inclined nature of the nebular shell, meaning that the 
slit cuts across a region where only the blue-shifted side of 
the nebula is present. However, projection effects, surface 
brightness and shell thickness variations could produce an 
apparently open velocity ellipse from a closed shell, due to 
the sensitivity limits of the observations. 

The intrinsic nebular structure and variation in bright- 
ness across the nebula are discussed further in Sections 13.21 
and [Til 



3.2 Spatio-kinematical modelling of Abell 41 

A spatio-kinematical model, corresponding to the simplest 
three-dimensional structure consistent with the large-scale 
nebular [Nil] 6584 A emission features, has been derived for 
A 41. [Nil] 6584 A emission was selected rather than Ho? 
emission due to its lower thermal broadening and its shell- 
like distribution (as discussed in section ITT]) . The modelling 
was performed in order to confirm the bipolar nature of 
the nebular shell and to constrain the inclination angle of 
this shell, for comparison with the inclination of the central 
binary (MT Ser). The model was developed using SHAPE 
ISteffen fc Lopej2006l ) and by assuming a Hubble- type flow, 
where expansion velocity is radial and proportional to the 
distance from the centre of the nebula. The model parame- 
ters (dimensions, shape, expansion scale velocity and incli- 
nation) were manually varied over a wide range of values 
and the results compared by eye to both spectral obser- 
vations and imaging, until a best-fit was found. This best 
fit model comprises a bipolar shell waisted by an equato- 
rial ring with an expansion velocity of ~40 km s~^. The 
model nebula is slightly asymmetric in that the northern 
lobe is shortened by 15%, has a narrower opening angle 
and has a slight shear with respect to its southern coun- 
terpart. No symmetric model could be found to reproduce 
the observed PV arrays. The nebular inclination angle, as 
defined by the un-sheared southern lobe, is determined to be 
66 ° ± 5° (in excellent agr eement with the value determined 
bv lPollacco fc Belli . Il997l . by deprojecting the nebular ring). 
The model nebula is shown at the observed orientation in 
Figure [S] and at an inclination of 90°, to highlight the asym- 
metry, in Figure|6l The synthetic PV arrays are shown, along 
with their observed counterparts, in Figure |3l 

It is not unheard of for bipolar nebulae to show asymme- 
try between the opposi ng lobes, both in the exten sion of the 
lobes (e.g. NGC 6881: iRamos-Larios et al.l 120081) and their 
opening angles (e.g. OH231.8-I-4.2: Isanchez Contreras et al.l 
I2OO4) . In comparison, the level of asymmetry shown by A 41 
is very low. One possible explanation for the asymmetry is 
discussed in Section [3.41 

Also of note, the equatorial expansion velocity, of 
~40 km s~^, is unusually large for a bipolar nebula, though 
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Figure 5. The synthetic SHAPE model for the Nil emission from 
A 41 displayed at the same scale as those images in Figure [T] 

not e xceptionally so (|Solf fc Weinbergeil 1 19841 . IWeinbergej 
1 19891 '). 

3.3 Systemic velocity and kinematical age 

Comparison of synthetic model spectra to their observed 
counterparts provides an unambiguous measure of the neb- 
ular systemic heliocentric radial velocity (Vsys), unaffected, 
for example, by brigh tness variations or nebular asymme- 
try (jjones et al.l 120101 ). Using the best-fit model described 
in Section [3.21 Vsys is determined to be 30 ± 5 km s^^ in 
good agreement with the value of 30 km s~^ determined by 
[Scaulicu ct al] l\199^ ). Similarly, the nebular expansion ve- 
locity, determined by the kinematical modelling, can be used 
to calculate a kinematical age for the nebula. This, however, 
requires the distance to the nebula to be known. The dis- 
tance to A 41 is a matter of so me debate with values in the 
literature rang ing from ~1 kpc (iGrauer fc Bondll 19831 ') up to 
9.0 ± 0.4 kpc (jShimanskii et ahboOSh . this probably resuhs 
from the notorious variation in results fr om different meth - 
ods of distance determination (see e.g., lGurzadvanlll997l ). 
Therefore, rather than favour one particular distance esti- 
mate over another we quote a kinematical age per kiloparsec 
of ~ 800 years kpc~^. 

3.4 Evidence of ISM interaction? 

The asymmetries in both brightness and shape (discussed 
in Sections 13.11 and 13. 2[) could be thought of as strong ev- 
idence for interaction with the interstellar medium (ISM). 
Consider, if the nebula were moving through the ISM with 
the northern lobe at the leading edge, it would be reason- 
able to expect to see this lobe brightened (through shock 
excitation) and less extended (due to the greater drag) with 
respect to its relatively unimpeded southern counterpart. 
This brightening and lesser extension is precisely what is 
borne out by the data and subsequent modelling. Similarly, 
if the motion of the nebula, with respect to the ISM, were 
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Figure 6. The synthetic SHAPE model for the Nil emission from 
A 41 displayed at an inclination of 90°. 



not along the symmetry axis but slightly offset, then this 
would also offer an explanation for the "shear" seen in the 
northern lobe (i.e. the shear is as a result of the eastern side 
of the lobe being impeded by the ISM more than the west- 
ern side) , and indeed the brightness contrast across that lobe 
(eastern edge brightened with respect to the western edge). 

We do not see any evidence of a bow-shock in any of 
our observations, which might be expe cted to be particularl y 
prevalent in the [Olll]5007 A image l|'Wareing et al.ll200'?l '). 
However, the relatively low levels of asymmetry in the neb- 
ula indicate that any interaction is fairly weak, implying a 
slow relative velocity consistent with no bow shock. 

The ratio of [Sii]6717-h6731 A to [Nii]6584 A provides 
a good tracer for shock excitation which might be expected 
to ar ise from interaction with the ISM (IPhillips fc Cuesta 
[1998'). However, the close proximity of the Moon to A 41 
during the [Sll]6717-(-6731 A imaging detailed in Section 
[2I and the uncertain nature of the subtraction of this 
background, means we are unable to accurately assess the 
[Sll]6717-f6731 A to [Nll]6584 A ratio of this nebula. Simi- 
larly, given the proper motion and Vsys of A 41 it would be 
possible to determine the motion of the nebula relative to 
th e local ISM (assumed to be due to Galactic rotation, as 
in lMeaburn et aLll2009l ). The direction of this motion could 
then be used to assess the validity of this hypothesis, unfor- 
tunately we are unable to perform this analysis as no proper 
motion measurements exist for this nebula. 

Also, of interest is the apparent location of the central 
star, which one would expect to be offset towards the leading 
edge of the nebula by any interaction between nebula and 
ISM (in the case of A 41 offset roughly towards the North 
of the nebula). In fact, in our imagery it is slightly offset to 
the south of the centre of the nebular ring, at odds with this 
hypothesis. 
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4 CONCLUSIONS 

Using high-resolution longsht spectroscopy and deep imag- 
ing, a spatio-kinematical model of A 41 has been developed 
which clearly shows that the nebula is aligned with the bi- 
nary central system exactly as predicted by current theo- 
ries of PN shaping by binary central stars. This is only the 
second neb ula to have this link observationally constrained 
(after A 63. lMitchell et al.ll2007l '). The kinematical data con- 
firm A 41 exhibits a waisted, bipolar structure, with some 
small deviations from perfect axisymmetry. The presence of 
an equatorial ring is also confirmed, adding further weight 
to the link between r ing structures and cen tral binary stars 
as commented on bv iMiszalski et al] (|2009l ). 

The data indicate that A 41 may be experiencing an 
interaction with the ISM and this too may be affecting its 
shape and brightness. 

Further kinematical investigations, such as the one pre- 
sented in this paper, coupled with in-depth studies of the 
CSPN of other PNe with confirmed close-binary central 
stars, are necessary to investigate the full extent of the in- 
fluence of central star binarity on PN nebular shaping. Only 
once a significant statistical sample has been acquired can 
generalisations be made about the role of CSPN binarity in 
PN evolution. 
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